Clinical significance and correlation of PD-L1, B7-H3, B7-H4, and TILs in pancreatic cancer

Background B7 molecules play significant roles in regulating tumor immunity, but their expression patterns and immuno-biological correlations in pancreatic cancer (PaCa) have not been fully discussed. Methods RNA-sequencing data of B7 molecules of PaCa samples in the Cancer Genome Atlas (TCGA) dataset was downloaded from the UCSC Xena to assess the expression, correlation, and mutation of the B7 family in PaCa. Next, two PaCa tissue microarrays (TMAs, Cat. HPanA150CS02 and HPanA120Su02) were obtained from Outdo BioTech (Shanghai, China). To detect the expression levels of PD-L1, B7-H3 and B7-H4, immunohistochemistry (IHC) staining was performed on these TMAs. Results Most B7 molecules, including B7–1, B7–2, PD-L1, B7-DC, B7-H2, and B7-H5 exhibited similar expression patterns, but B7-H3, B7-H4, B7-H6, and B7-H7 showed outlier expression patterns compared with other B7 molecules. Besides, B7 molecules were genetically stable and exhibited low alteration frequency. IHC staining indicated PD-L1, B7-H3, and B7-H4 were up-regulated in PaCa tissues and showed uncorrelated expression patterns. Furthermore, high expression of PD-L1 and B7-H3 indicated poor-differentiated grades in PaCa. PD-L1 was positively, but B7-H4 was negatively correlated with CD8+ TILs infiltration in PaCa. Moreover, combined PD-L1 and B7-H4 expression was a novel subtyping strategy in PaCa, namely patients with both high PD-L1 and B7-H4 expression exhibited decreased CD8+ TILs infiltration in tumor tissues. Conclusion Overall, we systemically analyzed the expression patterns of B7 molecules and proposed a novel subtyping strategy in PaCa. Patients with both high PD-L1 and B7-H4 expression exhibited the immuno-cold phenotype, which may be not suitable for immunotherapy. Supplementary Information The online version contains supplementary material available at 10.1186/s12885-022-09639-5.

worst prognosis and accounts for almost as many deaths as new cases [2]. In contrast to most other cancer types, improvement of the therapeutic options and prognosis over the last decades has been quietly limited [3,4].
Immune checkpoint blockade (ICB) is inspiring immunotherapy transforming the standard of treatment for several advanced tumors [5]. Representative immune checkpoint inhibitors are monoclonal antibodies that interfere with the interaction between PD-1 and PD-L1 inhibitory proteins expressed on the surface of T cells and tumor cells, respectively. Unfortunately, the response rate to ICB monotherapy remains low in PaCa [6,7]. A recent phase 2 clinical trial combining durvalumab plus tremelimumab exhibited a mere objective response (OR) rate of 3.1% [8]. Therefore, effective immune-related biomarkers are urgent to identify patients who will benefit from ICB and to stratify patients into optimal immunotherapy strategies.
Tumor mutation burden (TMB) and T cell inflamed gene expression profile are useful biomarkers for response to ICB in cancers [9]. In clinical practice, the detection of PD-L1 has been considered a predictive biomarker for anti-PD-1/PD-L1 therapies. Patients with high PD-L1 expression tended to be highly responsive to PD-1/PD-L1 blockade [10,11]. In addition to PD-L1, B7-H3 and B7-H4 are highly expressed in most tumor tissues, making them attractive candidate immunotherapeutic targets and biomarkers [12,13]. Furthermore, the combination of PD-L1 and B7-H4 expression is a novel classification strategy in glioma, namely the high B7-H4-expressed samples could be considered supercold gliomas with significantly deficient in tumor-infiltrating lymphocytes (TILs) [14]. However, the expression patterns of PD-L1, B7-H3, and B7-H4 in PaCa have not been fully discussed.
Given the evidence and the clinical need to explore immune-related biomarkers to stratify patients into optimal immunotherapy strategies, the current study aimed to systematically characterize the expression patterns and immuno-biological correlations of B7 family members, especially PD-L1, B7-H3, and B7-H4. Besides, based on the expression levels of PD-L1 and B7-H4, we proposed a novel subtyping strategy in PaCa, namely patients with both high PD-L1 and B7-H4 expression exhibited the immuno-cold phenotype, which may be not suitable for immunotherapy.

Acquisition of TCGA data
Normalized RNA-sequencing data of PaCa samples in the Cancer Genome Atlas (TCGA) dataset was obtained from the UCSC Xena (https:// xenab rowser. net/ datap ages/) or the cBioPortal portal (http:// www. cbiop ortal. org/ datas ets) [15]. Patients with missing or insufficient data were excluded from this research. The genetic alteration data of B7 molecules was downloaded from the cBi-oPortal portal.

Clinical specimens
Two PaCa tissue microarrays (TMAs, Cat. HPanA150CS02 and HPanA120Su02) were obtained from Outdo BioTech (Shanghai, China). The HPanA150CS02 TMA contained 78 PaCa tissues and 72 paired adjacent tissues, which was used as the exploring cohort for the detection of B7 molecules expression and the training cohort for the PD-L1/B7-H4-based stratified strategy. The HPanA120Su02 TMA contained 66 PaCa tissues and 54 paired adjacent tissues, which was used as the validated cohort for the PD-L1/B7-H4-based stratified strategy. Detailed clinico-pathological characteristics of the cohorts were also provided by Outdo BioTech (Shanghai, China). Ethical approval for the study of the TMAs was granted by the Clinical Research Ethics Committee, Outdo BioTech (Shanghai, China).

Statistical analysis
R 4.1.0, SPSS 26.0, and GraphPad Prism 6.0 were applied as the main tools for the statistical analysis and figures exhibition. Differences between the two groups were analyzed using Student's t-test or Mann-Whitney test for continuous variables and the chi-square test for categorical variables. Differences between multiple groups were analyzed using one-way ANOVA test or Kruskal-Wallis analysis for continuous variables. Correlation analysis was evaluated by Pearson correlation analysis. Log-rank test was conducted to evaluate the difference between the survival curves. For all statistical analyses, P ≤ 0.05 was deemed to be statistically significant.
We next validated the protein expression patterns of B7 molecules in PaCa using IHC staining. Considering the similar expression patterns of B7-1, B7-2, PD-L1, B7-DC, B7-H2, and B7-H5, we selected PD-L1 among these molecules for further validation. Besides, B7-H3 and B7-H4 were also selected for validation due to their relatively clear functions compared with B7-H6 and B7-H7. After analysis of IHC staining, we found that PD-L1, B7-H3, and B7-H4 proteins were all up-regulated in PaCa tissues compared with para-tumor tissues ( Fig. 2A-D). Moreover, similar to the results of the TCGA dataset, these three molecules were not correlated with each other (Fig. S1A-C). Overall, PD-L1, B7-H3, and B7-H4 were overexpressed in PaCa and showed uncorrelated expression patterns.

Association between expression levels of B7 molecules and clinic-pathological parameters in PaCa
Subsequently, we analyzed the association between expression levels of B7 molecules and clinic-pathological parameters. As shown in Table 1, PD-L1 and B7-H3 protein expression levels were not associated with age, gender, T stage, N stage, M stage, and clinical stage, but associated with tumor differentiation (Table 1). However, B7-H4 protein expression level was not associated with any clinic-pathological parameters in the cohort (Table 1). We next compared the expression levels of PD-L1, B7-H3, and B7-H4 in tumor tissues with various differentiated grades. The results showed that PD-L1 and B7-H3 proteins were up-regulated in poor-differentiated tissues, but B7-H4 protein expression was similar in welland poor-differentiated tissues (Fig. 3A-D). Moreover, the mRNA data from the TCGA dataset was also in accordance with the current results ( Fig. 3E-G). Due to the deficiency of survival data in the current cohort, we utilized the survival data from the TCGA dataset to assess the prognostic values of these three B7 molecules. However, the results showed that mRNA expression of single B7 molecules had limited prognostic value (Fig. S2A-C). Taken together, these results suggested that high expression of PD-L1 and B7-H3 indicated poor-differentiated grades in PaCa.

PD-L1/B7-H4 stratify infiltrating levels of TILs in PaCa
Given the correlation between PD-L1 and B7-H4 protein expression and TILs infiltration, we next tried to develop a novel subtyping strategy to estimate immunophenotype in PaCa. We divided PaCa patients into four subgroups, namely PD-L1 high -B7-H4 high , PD-L1 high -B7-H4 low , PD-L1 low -B7-H4 high , and PD-L1 low -B7-H4 low groups (Fig. 5A). According to previous research, high PD-L1 expression was correlated with increased CD8+ TILs infiltration and was the reliable biomarker for immuno-hot tumors [21]. However, the results showed that decreased CD8+ TILs infiltration in the PD-L1 high -B7-H4 high group compared with the PD-L1 high -B7-H4 low group (Fig. 5B). The HPanA120Su02 TMA contained 66 PaCa tissues was used as the validated cohort. In the HPanA120Su02 cohort, PD-L1 was positively correlated with CD8+ TILs infiltration but B7-H4 was not correlated with CD8+ TILs infiltration in PaCa ( Fig. S3A-C). The result of the HPanA120Su02 cohort was in accordance with the current result of the HPanA150CS02 cohort (Fig. 5C). Overall, these data suggested PaCa patients with both high PD-L1 and B7-H4 expression exhibited the immuno-cold phenotype, which may be not suitable for immunotherapy.

Discussion
Increasing numbers of clues suggest that the dys-regulated expression of B7 molecules results in the suppressive tumor microenvironment (TME) and immune escape [22]. The B7 family consists of 10 members, namely B7-1 (CD80), B7-2 (CD86), B7-H1 (CD274/ PD-L1), B7-DC (PDCD1LG2), B7-H2 (ICOSLG), B7-H3 (CD276), B7-H4 (VTCN1), B7-H5 (VSIR), B7-H6 (NCR3LG1), and B7-H7 (HHLA2) [23]. In this research, we first systemically explored the expression patterns of these B7 molecules in PaCa. B7-1, B7-2, PD-L1, B7-DC, B7-H2, and B7-H5 exhibited similar expression pattern, which was highly correlated with each other, but B7-H3,  B7-H4, B7-H6, and B7-H7 showed outlier expression patterns compared with other B7 molecules. Besides, B7 molecules were genetically stable with low alteration frequency in PaCa, suggesting that genetic alterations might not be critical threats to human cancers. In recent years, the dominating roles of the B7 molecules in regulating tumor immunity have been widely concerned, especially PD-L1, B7-H3, and B7-H4, which have become hot research objects [19,24]. To study their respective roles in the TME, the expression profiles and biological roles of PD-L1, B7-H3, and B7-H4 in cancers have been well-concluded. For example, Qiu et al. analyzed the expression and prognostic values of these B7 molecules and found that B7-H3 was a negative prognostic indicator in small cell lung cancer [25]. In this study, PD-L1, B7-H3, and B7-H4 were applied as research objects to be focused on, and their expression, clinical significance, and correlation with TILs in PaCa were explored. The critical results of the current study may In PaCa samples in the current cohort, we found that PD-L1, B7-H3, and B7-H4 were all overexpressed in PaCa tissues compared with para-tumor tissues, which was in accordance with other studies [26][27][28]. Besides, the protein levels of these three molecules were not correlated with each other as well. Moreover, PD-L1 and B7-H3 were up-regulated in poor-differentiated tissues, but B7-H4 expression was similar in well-and poor-differentiated tissues. However, survival analysis showed that single B7 molecules had limited prognostic value. Furthermore, we also checked the correlation between these B7 molecules expression and TILs levels in PaCa in this research. PD-L1 was positively, and B7-H4 was negatively correlated with CD8+ TILs infiltration, but B7-H3 was not correlated with CD8+ TILs infiltration in PaCa.
Growing numbers of research have reported that mutually-exclusive or co-expressed patterns of B7 molecules could predict inflamed or non-inflamed TME in multiple human cancers [14,29]. For example, B7-H4 is negatively correlated with PD-L1 and identifies immuno-cold tumors in glioma [14]. Considering the various correlations between PD-L1 and B7-H4 expression and TILs infiltration, we tried to develop a novel subtyping strategy to estimate immuno-phenotype in PaCa based on PD-L1 and B7-H4 expression. High PD-L1 expression has been considered a reliable biomarker for high response to immunotherapy due to its high correlation with effector TILs [21,[30][31][32]. First, PD-L1 was not correlated with B7-H4 expression in PaCa in the current cohort. Importantly, the results of our study indicated decreased CD8+ TILs infiltration in PaCa patients with both high PD-L1 and B7-H4 expression, and the result was also validated in another independent HPanA120Su02 cohort. Collectively, PaCa patients with both high PD-L1 and B7-H4 might be dismissed from the immunotherapy in clinical practice. The current research also had several shortcomings. First of all, the expression, correlation, and mutation of B7 family in PaCa were investigated using the public TCGA cohort, which were not validated using the in-house cohort. In addition, there may be a certain discrepancy between mRNA and protein level expression, and RNA-sequencing data can not distinguish the cell types. Furthermore, although we proposed a novel subtyping strategy in PaCa, it was not validated in any immunotherapy cohort in clinical practice.